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Una manera de impulsar proyectos con base en energías 
renovables, como los sistemas hidroeléctricos a pequeña 
escala, es determinando los lugares viables para el 
emplazamiento de dichos proyectos y la energía que podría 
explotarse. En este trabajo se expone un método de sencilla 
aplicación para estimar el potencial hidroenergético en 
cuencas hidrológicas. La bondad principal del método 
se debe al empleo de software libre para la modelación 
hidrológica continua. Se utiliza el software ArcGis 10.1 para 
realizar la delimitación y el procesamiento cartográfico de las 
subcuencas hidrológicas y el trazo de perfiles longitudinales 
de cauces. Los criterios establecidos para la selección de 
sitios disponibles fueron los siguientes: pendiente del río 
de al menos 3% y el desnivel topográfico de tres metros o 
más. Se elige un periodo común de precipitación y caudal 
medio diario para calibrar y validar el modelo hidrológico 
Moisture Accounting (SMA). El 
método SMA es una herramienta adecuada para conocer 


semidistribuido Soil 


el caudal medio en cuencas no aforadas. Se encontraron 
578 sitios con un potencial hidroenergético bruto de 
49.2 MW. La metodología empleada es práctica, porque 
requiere conocimientos básicos de sistemas de información 
geográfica y del manejo de HEC-HMS. La evaluación 
económica muestra inversiones iniciales de $3,672,636.00 
pesos mexicanos en promedio. Ante tales costos, puede ser 
factible invertir, si se considera una planeación a largo plazo. 
El Q,, es una buena alternativa para el diseño de pequeñas 
centrales hidroeléctricas al hilo de agua. 


Palabras clave: potencial hidroenergético, pequeñas 
centrales hidroeléctricas (PCHs), modelo Soil Moisture 
Accounting (SMA), curva de duración de caudales, río Alto 
Amacuzac. 


Abstract 


Meza-Prieto. B., E; Aparicio, J. (January-February, 2018). 
Evaluation of small-scale hydropower potential in watersheds 
using a rainfall-runoff model. Water Technology and Sciences 
(in Spanish), 9(1), pp. 69-87, DOI: 10.24850/¡-tyca-2018-01-05. 


One way of promoting renewable energy-based projects, such as the 
small-scale hydroelectric systems, is to determine viable placements 
for them, as well as the exploitable energy. This work presents a simple 
method to estimate the hydropower potential in watersheds. The 
main advantage of this method is that the continuous hydrological 
modeling is performed using free software. The delimitation and 
cartographic processing of the sub basins, as well as the longitudinal 
profile tracing of the riverbed is performed using the Arcgis 10.1 
software. The criteria for the available sites selection were established 
as: river slope of at least 3%, and vertical drop of 3 meters or more. 
A common precipitation and average daily flow period is chosen in 
order to calibrate and validate the Soil Moisture Accounting (SMA) 
semi-distributed hydrological model. The SMA method is a suitable 
tool to estimate the mean streamflow in ungauged watersheds. In 
the Alto Amacuzac River Basin, 578 sites with a gross hydropower 
potential of 49.2 MW were found. The methodology used is practical 
because it only needs basic knowledge of geographic information 
systems and uses HEC-HMS. The economical evaluation shows 
initial investments of $3,672,636.00 Mexican pesos on average. 
With these costs, investments may be feasible as long as a long-term 
plan is considered. The streamflow Q,, is a good alternative to design 
run-of-river small hydropower plants. 


Keywords: hydropower potential, small hydropower plants, Soil 
Moisture Accounting (SMA) model, flow-duration curve, Alto 
Amacuzac River. 
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Introducción 


La generación de energía hidroeléctrica a peque- 
ña escala es una actividad de gran importancia 
en el mundo y se ha practicado desde hace más 
de doscientos años (Fritz, 1984). Además de 
ser ecológicamente más conveniente (Castro, 
2006; APER, 2007), las ventajas económicas de 
las hidroeléctricas sobre los sistemas de energía 
térmica son muy considerables, pues el costo de 
generación es menor comparado tanto con los 
sistemas de generación eléctrica nuclear como 
con los que funcionan a partir de combustibles 
fósiles, como diésel y gas natural (Calderón, 
2017). 

Una alternativa para las plantas de gene- 
ración hidroeléctrica a pequeña escala son las 
construidas al hilo del agua, pues su impacto 
ambiental es despreciable, comparado con las 
grandes centrales, las cuales se caracterizan por 
generar energía de base a todas horas, aunque 
tienen una carga constante con limitada adap- 
tación a la demanda de energía y no permiten 
regular las inundaciones ni la variación del flujo 
en el río (Ortiz, 2011). 

El potencial eléctrico de pequeñas y micro 
centrales hidroeléctricas ha sido evaluado en 
países como Estados Unidos y Brasil. En Esta- 
dos Unidos se tiene un potencial de energía total 
bruto anual en este tipo de centrales de cerca de 
171 055 MW, pero sólo están disponibles 101 341 
MW. En el caso de Brasil, se estimó un potencial 
para pequeñas centrales hidroeléctricas de 237 
870 MW, de los cuales 4 822 se encuentran en 
sitios ya desarrollados y 50 231 están en áreas 
de exclusión, y por tanto sólo pueden utilizarse 
182 817 MW (Hall, 2011). 

De acuerdo con Alemán et al. (2014), en Mé- 
xico existen algunas barreras para el desarrollo 
de energías renovables, entre las que se pueden 
mencionar: la estrategia energética, que se basa 
en la evaluación de los recursos de energía 
renovable a corto plazo y no a largo plazo; la 
tecnología, por la insuficiente inversión en la 
exploración de los recursos de energía reno- 
vable que no han sido considerados, como los 
de energía geotérmica y la undimotriz, y en el 
incremento de la capacidad de transmisión en 
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áreas con alto potencial de recursos de energía 
renovable; y los incentivos, pues si bien la Comi- 
sión Reguladora de Energía (CRE), la Comisión 
Federal de Electricidad (CFE) y la Secretaría 
de Energía (Sener) de México los ofrecen para 
estimular la participación del sector privado en 
la generación de energía de recursos renovables, 
ello no es suficiente para el desarrollo de todos 
los esquemas de generación eléctrica. 

Por otro lado, la generación por medio de 
pequeñas hidroeléctricas enfrenta otras barre- 
ras, como la falta de evaluación confiable del 
potencial de generación y la incertidumbre 
debida a la escasez de información meteoroló- 
gica e hidrométrica adecuada. Aunado a esto, 
existe otra limitante en el aprovechamiento de 
agua para generar electricidad debido a que 
México tiene grandes zonas áridas y esto lo 
hace vulnerable a las sequías, según la Comi- 
sión Nacional del Agua (Conagua, 2011) por las 
razones siguientes: a) hay periodos de tiempo 
en los que disminuye la producción de energía 
hidroeléctrica por falta de agua; b) el régimen de 
lluvias sufre alteraciones tanto espaciales como 
temporales; c) la evapotranspiración se incre- 
menta en los meses de marzo a junio, y d) se 
deben satisfacer las demandas de agua potable 
y riego antes que las de energía. 

En el caso de generación a gran escala, los 
sitios aprovechables son limitados en número 
y éstos no necesariamente serán rentables a 
largo plazo debido al cambio climático, que 
producirá una variación temporal en el régimen 
de lluvias en verano; en efecto, la precipitación 
podría disminuir en toda la república hasta en 
6% en primavera-verano en el periodo 2015- 
2039 de acuerdo con el escenario RCP4.5 del 
Panel Intergubernamental de Cambio Climático 
(IMTA, 2015). Por ello, para usar la energía hi- 
droeléctrica en el futuro, es conveniente realizar 
una planeación de los recursos a largo plazo, 
considerando, entre otros aspectos, la genera- 
ción en pequeña escala, es decir, construyendo 
pequeñas centrales hidroeléctricas (PCHs). 

El objetivo general de este trabajo es propo- 
ner un método para la estimación del potencial 
hidroenergético en cuencas hidrológicas y 
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localizar los sitios factibles para construir mi- 
cro, mini y pequeñas centrales hidroeléctricas. 
Sus objetivos particulares son: 1) implementar, 
calibrar y evaluar el modelo Soil Moisture Ac- 
counting (SMA) (Bennett ¿e Peters, 2000) como 
herramienta para obtener el caudal medio dis- 
ponible en cuencas hidrológicas para generar 
energía eléctrica; 2) definir lugares apropiados 
en las subcuencas del río Alto Amacuzac para 
construir centrales hidroeléctricas al hilo de 
agua y calcular su potencial; y 3) realizar la 
evaluación de los costos aproximados de insta- 
lación de una PCH y revisar su factibilidad para 
un caso específico. 


Metodología 


Como parte de la metodología se realiza un 
análisis cartográfico de las subcuencas hidro- 
lógicas mediante el software ArcGIS 10.1 (ESRI, 
2017). Por otro lado se analiza la información 
climatológica e hidrométrica, y posteriormente 
se calibra y valida el modelo hidrológico SMA 
usando el software HEC-HMS 4. El programa 
HEC-HMS tiene varias alternativas para la 
simulación de eventos de lluvia de corta du- 
ración y simulación continua. Para determinar 
el volumen total de agua disponible para 


escurrimiento superficial, el programa cuenta 
con el modelo Soil Moisture Accounting (SMA) en 
el caso de simulación continua (USACE, 2017). 
En la figura 1 se muestra el diagrama que ilustra 
el procedimiento propuesto para evaluar el po- 
tencial hidroenergético en cuencas hidrológicas. 


Análisis cartográfico 


El análisis cartográfico se realizó a partir del 
modelo digital de terreno de escala 1:50 000, 
del conjunto de datos vectoriales de uso de 
suelo y vegetación, escala 1:250 000 serie IV 
(INEGI, 2010) y del conjunto nacional de datos 
vectorial edafológico, escala 1:250 000 serie II 
(INEGL 2007), obtenidos del Instituto Nacional 
de Estadística y Geografía (INEGI). Dicho aná- 
lisis consistió en determinar las características 
fisiográficas de las subcuencas y los perfiles 
longitudinales de las corrientes perennes. Se 
identificaron los tipos de vegetación y suelo, y 
los elementos de exclusión. 


Determinación de las características fisiográficas de 
las subcuencas 


El río Amacuzac es un afluente del río Balsas 
y se origina en las faldas del volcán Nevado 
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Figura 1. Procedimiento para estimar el potencial disponible y costos. 
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de Toluca, a una altitud aproximada de 2 600 
msnm dentro de San Martín Tequesquipan, Es- 
tado de México, en donde se le conoce como río 
Texcaltitlán. La longitud total del río Amacuzac 
hasta la confluencia con el río Balsas es de 253 
km aproximadamente (IMTA, 2012). El área de 
estudio se encuentra ubicada en la cuenca río 
Alto Amacuzac, cuya superficie es de 2 611.98 
km, dentro de la cual se encuentran los ríos 
principales Texcaltitlán, Chontalcuatlán y San 
Jerónimo (INEGI, 2017). 

Se definieron los sitios propicios para cons- 
truir pequeñas centrales hidroeléctricas en esta 
cuenca considerando los siguientes criterios: 1) 
la pendiente del tramo de río debe ser por lo 
menos de 0.03 (3 %) debido a la resolución del 
modelo digital de terreno; y 2) un desnivel mí- 
nimo por tramo de tres metros debido a que las 
turbinas comerciales funcionan con más de dos 
metros de carga (Ortiz, 2011). Posteriormente se 
delimitaron las subcuencas y se obtuvieron sus 
características principales, como área, longitud 
del cauce principal, pendiente media del cauce 
principal y tiempo de concentración necesarias 
para el uso del modelo SMA. 


Determinación de perfiles longitudinales 


Se trazaron los perfiles longitudinales de las 
corrientes tributarias que coincidían con las co- 
rrientes perennes de las subcuencas hidrológicas 
(INEGL, 2017), con el fin de determinar el des- 
nivel topográfico y la distancia horizontal del 
cauce, y luego calcular una pendiente media. En 
la figura 2 se muestra un perfil tipo obtenido con 
la herramienta 3D-Analyst. Para ello se convier- 
ten los perfiles en 2D a 3D utilizando la rutina 
siguiente: ArcToolbox-3D Analyst Tools-Functional 
Surface-Interpole Shape. 


Identificación de los tipos de vegetación y suelo 


Se determinaron el tipo de vegetación y tipo 
de suelo en el área de estudio para definir el 
valor de los parámetros iniciales que dependen 
de las características geomorfológicas de las 
cuencas, tomando como base los reportados 
en otros trabajos. En la figura 3 se muestra el 
tipo de vegetación característica y uso de suelo 
de la cuenca Alto Amacuzac. El tipo de suelo 
dominante de acuerdo con la clasificación 
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Figura 2. Perfil longitudinal del cauce principal del río Amacuzac. 
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Tipo de vegetación 
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Figura 3. Tipo de vegetación y uso de suelo. Fuente: elaboración propia con datos de INEGL. 


del Servicio de Conservación de Recursos 2004). En la figura 4 se presenta la ubicación 
Naturales de Estados Unidos —antes Soil Con- de los tipos de suelos principales en la cuenca 
servation Service (SCS) (Aparicio, 2010) — se Alto Amacuzac. La clase NA corresponde a 
obtuvo a partir de la información edafológica las áreas donde están ubicadas las manchas 
digital con clases de textura 1, 2 y 3 (INEGI, urbanas. 
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Tipo de suelo 
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Figura 4. Tipo de suelo y área urbana. Fuente: elaboración propia con datos de INEGL. 


Ubicación de los elementos de exclusión 


Como elementos de exclusión se consideran las 
superficies en las que no se puede aprovechar 
el agua para generación hidroeléctrica, por 
ejemplo las áreas protegidas, proyectos de 
Obras existentes, zonas con baja precipitación 
y cuando el uso de suelo no lo permite, como 
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en el caso de áreas urbanas o agrícolas. En la 
figura 5 se muestra el área urbana, infraestruc- 
tura existente y las áreas naturales protegidas 
pertenecientes a los parques nacionales Nevado 
de Toluca y Desierto de Carmen o de Nixcongo, 
así como el Parque Estatal Santuario del Agua 
Temascaltepec, ubicados dentro de la cuenca 
Alto Amacuzac. 
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Figura 5. Ubicación de las áreas naturales protegidas y área urbana. Fuente: elaboración propia con datos de INEGI. 


Análisis de información climatológica e 
hidrométrica 


Para realizar el análisis de la información cli- 
matológica se utilizaron los datos registrados 
en la base de datos climatológica nacional 
CLICOM del Servicio Meteorológico Nacional, 
con información de precipitación hasta el año 


2012 (Conagua, 2017). La información hidro- 
métrica se obtuvo del Banco Nacional de Datos 
de Aguas Superficiales (Bandas) de la Red Hi- 
drométrica Nacional, actualizado por Instituto 
Mexicano de Tecnología del Agua (IMTA, 2012). 
Asimismo, se obtuvo la evapotranspiración 
mensual con el método Blaney-Criddle (Apari- 
cio, 2010), ya que éste utiliza datos disponibles, 
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como temperatura media mensual y horas de 
sol mensual. Se definió un periodo común de 
al menos 12 años de datos registrados de preci- 
pitación y caudal medio diario para usar en la 
calibración y validación del modelo SMA, pues 
se consideró el criterio tomado en la calibración 
y validación del modelo SMA en HEC-HMS, en 
el que se calibró el modelo con un periodo de 
datos de caudal de cuatro años, menor al usado 
en la validación de siete años. 


Revisión de información climatológica 


Se trazaron los polígonos de Thiessen (Apa- 
ricio, 2010) mediante el software ArcMap para 
obtener el área de influencia fija de cada es- 
tación climatológica necesaria para asignar 
un peso espacial a la precipitación registrada 
en el modelo meteorológico del software 
HEC-HMS. El uso de polígonos de Thiessen 
permite determinar áreas de influencia de las 
estaciones climatológicas para varias cuencas 
a la vez. Para esto se utilizaron 23 estaciones 
climatológicas, de las cuales cinco están ubi- 
cadas en Morelos, 14 en el Estado de México y 
cuatro en Guerrero. Se revisaron los datos de 
precipitación diaria y se dedujeron los datos 
faltantes para completar las series de datos 
históricos diarios de acuerdo con Aparicio 
(2010). Se encontraron más datos faltantes en 
el periodo 1991-2006 utilizados para realizar 
la validación de los parámetros en la cuenca 
Chontalcuatlán, pues de las nueve estaciones 
climatológicas, cuatro presentaron gran porosi- 
dad; por ejemplo, en Meyuca faltaban 12 meses 
de información diaria; en San Fco. Oxtotilpan, 
20 meses; en Coatepec Harinas, 33 meses y en 
Nevado de Toluca faltaban nueve meses. En la 
estación Troja Vieja se determinaron los datos 
faltantes de dos meses; en Coatepequito y Dos 
Bocas se rellenaron dos meses de precipitación 
diaria, y en El Mirador un mes. En la mayoría 
de los casos se estimaron los datos faltantes de 
los meses menos lluviosos y en todos los casos 
se observó una correlación lógica en los datos 
de precipitación mensual. 
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Análisis de caudales medios diarios 


Para hacer el análisis se emplearon cinco esta- 
ciones hidrométricas ubicadas en la cuenca Alto 
Amacuzac, sobre los cauces de los ríos Amacu- 
zac, San Jerónimo, Chontalcuatlán y Texcaltit- 
lán. No existen reglas generales para seleccionar 
los periodos de calibración y validación, pero 
algunos autores (e.g., Bottcher, Whiteley, James, 
£z Hiscock, 2012) recomiendan que el periodo 
de validación sea al menos del doble que el de 
calibración. En los casos de las cuencas Texcal- 
titlán y Chontalcuatlán se usaron cuatro años de 
datos históricos de caudales diarios, así como 8 
y 16 años, respectivamente, para la validación; 
los parámetros para la cuenca Coatepequito se 
calibraron con cinco años y se validaron con 
20 años de registros, y los parámetros para las 
cuencas Dos Bocas y Amacuzac se calibraron 
con seis años y se validaron con 16 años de 
registros de datos diarios. Esto se hizo con el fin 
de utilizar los datos de precipitación empleados 
en la validación de los parámetros, para hacer 
la simulación en las subcuencas hidrológicas, 
pues, como Ortiz (2011) menciona, las curvas de 
duración de caudales y de frecuencias se deben 
obtener con la información histórica existente 
superior a 10 años preferiblemente en el caso 
de PCHs. 

En la figura 6 se muestra la ubicación de 
las estaciones climatológicas e hidrométricas 
utilizadas para el análisis y los polígonos de 
Thiessen trazados. El nombre de la cuenca hi- 
drológica se asignó de acuerdo con el nombre 
de la estación hidrométrica de aforo. 


Calibración y validación del modelo 
hidrológico 


El modelo hidrológico está constituido por tres 
etapas de cálculo. En la primera se determina 
el flujo del agua de precipitación con el modelo 
SMA, ya sea que se evapotranspire, se infiltre, 
percole o escurra. Luego se utiliza el modelo 
de Clark para representar la transformación 
de la precipitación en exceso en escurrimiento 
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Figura 6. Ubicación de las estaciones climatológicas e hidrométricas. Fuente: elaboración propia con datos de INEGI y Conagua. 


superficial; finalmente se determina la contri- 
bución del flujo base en el hidrograma resul- 
tante mediante el modelo de depósitos lineales 
(Bennett, 1998). 

En el proceso de calibración es necesario 
definir valores adecuados para cada paráme- 
tro; entre mayor sea el número de parámetros 


involucrados en cada etapa de cálculo, mayor es 
el grado de dificultad para obtenerlos. En este 
caso, para cada cuenca hidrológica se deben 
establecer 25 parámetros: 17 pertenecientes al 
modelo SMA, dos pertenecientes al modelo 
de Clark y seis pertenecientes al modelo de 
depósitos lineales. 
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Se puede recurrir a varias estrategias para 
calibrar un modelo con varios parámetros. Si los 
parámetros tienen sentido físico, se establecen 
relaciones entre sus valores y las características 
geomorfológicas de las cuencas, utilizando 
expresiones matemáticas, tablas y rangos de va- 
lores orientativos propuestos por diferentes au- 
tores; si los parámetros no tienen sentido físico, 
pero dominan en ellos algunas condiciones de 
tipo matemático, en algunos casos hay valores 
que se pueden estimar por criterios subjetivos, 
comparándolos con los adoptados en estudios 
similares (Grupo de Emisarios Submarinos e 
Hidráulica Ambiental de la Universidad de 
Cantabria, 2004). En este estudio se obtuvieron 
los parámetros a partir de las características 
geomorfológicas, ya que varios parámetros 
dependen del tipo de vegetación; la pendiente 
del terreno; el uso y el tipo de suelo, como es 
el caso de la capacidad de almacenamiento del 
dosel forestal y la capacidad de almacenamiento 
superficial en depresiones. De igual manera se 
determinaron los parámetros relacionados con 
las características del suelo (textura, porosidad 
y permeabilidad), como la máxima capacidad 
de infiltración de agua en el suelo, la máxima 
capacidad de almacenamiento de agua en el 
suelo, la capacidad de almacenamiento en la zo- 
na de tensión del suelo y la máxima capacidad 
de percolación. 


Análisis de sensibilidad 


Para realizar el análisis de sensibilidad se pro- 
pusieron valores de los parámetros tomando 
como base los reportados en los trabajos de Ben- 
nett (1998) y el Grupo de Emisarios Submarinos 
e Hidráulica Ambiental de la Universidad de 
Cantabria (2004), y se realizó la simulación en la 
cuenca Texcaltitlán. Después se optó por multi- 
plicar por factores de 10 y de 0.2, cambiando uno 
a la vez los valores de los parámetros necesarios 
para realizar la simulación en el modelo SMA, 
excepto cinco valores que representan el alma- 
cenamiento en la vegetación, en superficie, en el 
suelo y subsuelo al inicio de la simulación, pues 
estos valores no afectan los resultados obtenidos 
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cuando se hace la simulación continua de varios 
años (Bennett, 1998). 


Calibración del modelo 


Para realizar la calibración se optimizaron los 
parámetros, utilizando las funciones objetivo 
“suma de errores absolutos”, “suma de resi- 
duos cuadrados”, “error porcentual en pico”, 
“error cuadrático medio de la raíz del pico- 
ponderado” y “error porcentual en volumen” 
(USACE, 2017). Debido a que en este caso inte- 
resa conocer el caudal de escurrimiento diario 
en las subcuencas, se consideró como criterio 
de aceptación un error menor de 10% entre los 
volúmenes observado y simulado. 


Validación del modelo 


Para determinar la bondad de ajuste del modelo 
hidrológico se usa el coeficiente de eficiencia de 
Nash-Sutcliffe (Nash 8 Sutcliffe, 1970). Según 
Moriasi et al. (2007), los valores del coeficiente 
de eficiencia de Nash-Sutcliffe entre 0 y 1 son 
aceptables, y el óptimo es el valor 1, mientras 
que un valor menor que 0 indica que la media 
de los valores observados es un mejor predictor 
que el modelo. 


Estimación del potencial disponible 


Para cada sitio se determinó la curva de dura- 
ción de caudales utilizando los datos obtenidos 
de caudal medio con el modelo calibrado en 
HEC-HMS. De la curva de duración de cauda- 
les se obtuvo el caudal mínimo que se puede 
asegurar con alta probabilidad, normalmente 
entre 85 y 95%. En este trabajo se utiliza el Q,, 
es decir, el caudal que se encuentra en el río el 
90% del año (UTP, 2010; Ayros € Salazar, 2011). 
Posteriormente se calculó la potencia hidráulica 
bruta utilizando la ecuación: 


P=9.81Q0H (1) 


donde P es la potencia bruta en kW, Q es el 
caudal en m*/s y H es la carga de presión debida 
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a la energía potencial en m, que en este caso es 
el desnivel topográfico (Ortiz, 2011). 


Análisis de costos 


El costo del kilowatt instalado es variable aun 
en centrales diseñadas con la misma potencia, 
pues depende de las condiciones propias del 
sitio y de las características técnicas del proyec- 
to. Con el fin de conocer el costo aproximado 
de la inversión que se realizaría al llevar los 
sitios disponibles a las siguientes etapas de 
planeación y construcción de las PCHs, se hizo 
la evaluación de los sitios disponibles utilizando 
información del costo por kilowatt instalado 
(en euros) reportado en Ortiz (2011), aunque en 
algunos países puede variar el costo debido a 
factores como el transporte y la importación de 
los equipos. El costo de instalación de la PCH 
depende naturalmente de la caída y la potencia, 
y las centrales con potencia menor de 250 kW 
son más costosas por unidad de potencia, como 
se puede observar en la figura 7. 


Se estimó el costo de cada PCH y la energía 
generada anualmente, considerando la eficiencia 
del grupo turbina-generador de 80%, pues se 
consideran turbinas hidráulicas tipo Michell 
Banki (Ortiz, 2011) de admisión parcial, por ser 
de bajo costo y ser adecuadas para desniveles 
bajos, además de ser muy estables de 1/3 a 3/3 
de carga. 

Los egresos anuales por los costos de opera- 
ción anual, reparación y mantenimiento anual 
de la planta e insumos para administración 
anual fueron considerados como el 3.5% del 
costo por kilowatt instalado de acuerdo con la 
distribución porcentual de los costos de cons- 
trucción de una PCH (Ortiz, 2011). El ingreso 
anual es por concepto de la venta de energía a 
un costo promedio de $1.45 por kWh, obtenido 
con los costos del consumo básico, intermedio 
y excedente de las tarifas domésticas 1, 1A, 1B, 
1C, 1D, 1E y 1F en temporada de verano y fuera 
de verano, que varían en el rango de $1.26 a 
$1.52, según la Comisión Federal de Electricidad 
de México (CFE, 2017). Con esta información se 
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Figura 7. Costo por kilowatt instalado para PCHs (Ortiz, 2011). 
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obtuvo el ingreso neto anual como la diferencia 
entre los egresos e ingresos anuales, información 
necesaria para analizar la rentabilidad de cada 
proyecto. 


Factibilidad de una PCH al hilo de agua 


Al realizar el estudio de factibilidad es difícil 
decidir cuál debe ser el caudal de diseño para 
las pequeñas centrales hidroeléctricas. Por ello, 
como ejemplo, se utiliza la curva de duración 
de caudales para determinar la viabilidad del 
proyecto con base en cada uno de los caudales 
presentes en el cauce durante el año en el sitio 
351, con código 7 749, ubicado en la cuenca 
Chontalcuatlán. La potencia de diseño es la 
potencia obtenida con la carga neta, es decir, la 
carga bruta o disponible menos las pérdidas en 
la tubería. La potencia mínima es la potencia a la 
cual funciona la turbina de manera aceptable; en 
este caso, dicha potencia se obtiene con el 16% 
del gasto de diseño, de acuerdo con la curva 
de rendimiento que se mantiene casi horizontal 
entre el 16 y 100% del caudal máximo de diseño 
(Osserberg, 2017). 

Se propone considerar acero inoxidable con 
un coeficiente de rugosidad n = 0.014 para deter- 
minar las pérdidas en la tubería, pues se buscó 
el diámetro de tubería con el que se perdiera 
menos de 10% de la carga disponible dada por 
el desnivel topográfico. 

Se presentan dos criterios para determinar 
la rentabilidad de los proyectos: la tasa interna 
de retorno (TIR) y el valor presente neto (VPN). 
La primera es la tasa efectiva anual con la que 
el valor presente neto de los costos es igual al 
valor presente neto de los beneficios (Fernán- 
dez, 2007). Para este caso se considera que la 
vida útil de cada proyecto es de 25 años, pues 
aunque la obra civil y la tubería forzada pueden 
tener una vida útil de unos 50 años, la vida útil 
de la turbina es por lo general de 35 años, y la 
del generador de 14 a 25 años, dependiendo de 
la potencia, por lo que el proyecto requeriría 
una nueva inversión en el generador o turbina 
si se sigue utilizando la pequeña central hasta 
los 50 años. 
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Resultados y discusión 


Calibración y validación del modelo 
hidrológico 


Análisis de sensibilidad 


Del análisis de sensibilidad se definió que la 
capacidad de almacenamiento del suelo (soil 
capacity); la capacidad del almacenamiento en la 
zona de tensión del suelo (tension zone capacity); 
la capacidad de almacenamiento del depósito 
subterráneo +1 (groundwater 1 storage coefficient); 
la capacidad de almacenamiento del depósito 
subterráneo +2 (groundwater 2 storage coefficient); 
la máxima capacidad de percolación del depósi- 
to subterráneo 42 (groundwater 2 percolation rate), 
y los coeficientes de almacenamiento del depó- 
sito subterráneo + 1 y 4 2 (groundwater 1 storage 
coefficient y groundwater 2 storage coefficient) son 
los parámetros que más influyen en la variabili- 
dad del caudal obtenido en la simulación, por lo 
que se debe poner especial atención al proponer 
los valores iniciales de estos parámetros para 
calibrar del modelo. 


Calibración del modelo 


Al realizar la simulación con las diferentes fun- 
ciones objetivo se observó mayor similitud entre 
los hidrogramas medido y simulado con la fun- 
ción “error porcentual en volumen”. También 
se advirtió que el modelo tiende a establecerse 
con el primer evento de lluvia y puede modelar 
con cierta precisión desde un año de datos con- 
tinuos, por lo que se puede deducir que no hay 
gran diferencia si se utiliza un periodo menor 
de datos en la calibración o en la validación. Los 
parámetros optimizados obtenidos en el proceso 
de calibración se muestran en el cuadro 1. 


Validación del modelo 


Al determinar la eficiencia del modelo para 
simular el volumen acumulado bajo el hidro- 
grama se encontraron muy buenos resultados, 
pues el coeficiente de Nash-Sutcliffe fue de 0.9 o 
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Cuadro 1. Parámetros optimizados para cada una de las cuencas hidrológicas. 


Núm. Parámetro Unidades | Texcaltitlan | Coatepequito | Chontalcuatlán | Dos Bocas | Amacuzac 
1 [Canopy Capacity mm 1.609 1.174 1.063 1.155 1.154 
2  |Canopy Initial Storage Percentage % 0.100 0.228 0.147 0.136 0.175 
3  |Clark Storage Coefficient h 846.090 780.060 780.100 780.020 780.070 
4 |Clark Time of Concentration h 2.110 5.980 8.210 10.100 11.130 
5  |Groundwater 1 Capacity mm 20.121 45.073 54.988 55.124 55.162 
6 Peal iS % 0.102 0.254 0.205 0.418 0.188 
7  |Groundwater 1 Percolation Rate mm/h 0.734 0.631 0.907 1.106 0.873 
8  |Groundwater 1 Storage Coefficient h 3 300.100 6 000.100 7 000.100 7.000.000 | 7.000.100 
9  |Groundwater 2 Capacity mm 16.651 38.063 50.138 50.000 50.167 
10 o AS % 0.102 0.253 0.197 0.084 0.183 
11 |Groundwater 2 Percolation Rate mm/h 2.118 1.566 0.918 0.894 0.721 

12  |Groundwater 2 Storage Coefficient h 3 522.100 4 750.100 5 999.700 6.000.000 | 6.000.200 
13  |Linear Reservoir GW 1 Coefficient h 3 300.100 6 000.100 7 000.000 6999.900 | 7.000.000 
14  |Linear Reservoir GW 1 Steps 2 15 18 28 30 
15  |Linear Reservoir GW 2 Coefficient h 3 522.100 4 750.100 6 000.100 6 000.000 | 6000.200 
16  |Linear Reservoir GW 2 Steps 2 15 18 28 30 
17 [Soil Capacity mm 2.613 2.874 2.090 2.184 2.175 
18 [Soil Infiltration Rate mm/h 2.266 1.463 0.769 1.186 0.839 
19  |Soil Initial Storage Percentage % 0.089 0.248 0.151 0.173 0.190 
20  |Soil Percolation Rate mm/h 0.531 0.445 1.111 0.755 0.689 
21 ¡Surface Capacity mm 1.127 1.152 1.075 1.190 1.153 
22  |Surface Initial Storage Percentage % 0.096 0.230 0.144 0.122 0.162 
23  |Tension Zone Capacity mm 0.541 0.749 0.574 0.606 0.647 


mayor en todas las cuencas hidrológicas, como diarios que los gastos pico. En la figura 8 se 


se puede observar en el cuadro 2, mientras que muestran los resultados de la validación de los 
para el caudal pico simulado, el modelo no parámetros para la cuenca Amacuzac. 


es tan eficiente, pues el caudal pico simulado 


representa menos de 50% del caudal pico ob- 
servado en las cinco cuencas. Sin embargo, el 
modelo es eficiente al representar los caudales 
medios diarios en un periodo de simulación de 
ocho años para la cuenca Texcaltitlán, 20 años 
para Coatepequito y 16 años para las cuencas 
Chontalcuatlán, Dos Bocas y Amacuzac, pues 
se observan valores entre “satisfactorio” y 
“adecuado” (Moriasi et al., 2007). 

Debido al carácter de este trabajo, es más 
importante representar adecuadamente los 
volúmenes escurridos y los caudales medios 


Estimación del potencial disponible 


En las 578 curvas de duración de caudales, los 
valores de Q,, variaron entre 0.020 y 3.76 m*/s 
para obtener el potencial bruto disponible. En la 
figura 9 se muestran curvas tipo seleccionadas al 
azar de las cuencas Texcaltitlán, Coatepequito, 
Chontalcuatlán y Dos Bocas. 

Se encontró que el potencial hidroenergéti- 
co total de la cuenca Alto Amacuzac es de 49 
227 kW distribuido en 578 sitios, de los cuales 
575 tienen un potencial menor de 1 MW. En 
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Cuadro 2. Eficiencia del modelo para simular el volumen acumulado, caudal pico y el caudal medio diario. 


Volumen acumulado Caudal pico Caudal medio diario 
Nash Caudal pico Caudal pico Nash- 
Cuenca , Ajuste observado simulado Sutcliffe Ajuste 
Sutcliffe 
(m/s) (m/s) 

Texcaltitlán 0.92 Muy bueno 26.3 7.5 0.610 Adecuado 
Coatepequito 0.95 Muy bueno 101.4 42.1 0.614 Adecuado 
Chontalcuatlán 0.89 Muy bueno 154.2 70.3 0.397 Satisfactorio 
Dos Bocas 0.98 Muy bueno 236.2 101.9 0.562 Adecuado 

Amacuzac 0.96 Muy bueno 506.6 147.3 0.518 Satisfactorio 
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Figura 8. Hidrograma obtenido para la cuenca Amacuzac, periodo de validación. 


el cuadro 3 se muestran las coordenadas de 
localización de algunos sitios disponibles; el 
desnivel o carga bruta, y la distancia horizontal 
en el tramo considerado; la pendiente media del 
tramo; el área de aportación al sitio elegido; el 
caudal presente en el cauce el 90% de año, y la 
potencia bruta disponible en los sitios elegidos. 
Dentro de los sitios localizados, el que tiene 
el potencial máximo es de 1 831.07 kW y el de 
menor potencial es de 1.43 kW. 

En la figura 10 se presenta un mapa con 
la distribución de los sitios, resaltando el po- 
tencial disponible por zonas de acuerdo con 
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la ubicación de los ríos y arroyos respecto al 
área urbana: El Chiquihuitero hasta el río de 
las Flores, Chiltepec y Palo Seco al norte de 
Texcaltitlán; Los Capulines en las localidades 
de Las Vueltas y Pachuquilla; Yerbabuena y 
El Cristo al oeste del poblado La Unión Riva 
Palacios; arroyo Florido (segunda sección), 
Ayatuxco, Arroyo Santiago y río El Alacrán 
cerca de Zacualpan; Acevedocla y San José 
cerca de la localidad Poder de Dios y el mu- 
nicipio de Tetipac; arroyo Los Tizantes, arroyo 
San Gaspar y río San Mateo al Norte de To- 
tolmajac; arroyo Los Cuervos, arroyo Tintojo 


Meza-Prieto £ Aparicio, Evaluación del potencial hidroenergético a pequeña escala en cuencas hidrológicas con un modelo lluvia-escurrimiento 


Cuenca Texcaltitlán 
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Figura 9. Curvas de duración de caudales tipo en las cuencas hidrológicas. 


Cuadro 3. Características de algunos de los sitios localizados. 


Sitio | Código | Desnivel | Distancia | Pendiente Esad a Potencia Coordenada | Coordenada 
(m) (m) Media pra A ruta (KW) | UTMY UTMX 
1 398 13.40 235.25 0.057 673.669 2.920 383.716 2072913.4713 | 439454.0853 
70 1554 68.68 957.68 0.072 160.969 0.900 606.385 2095089.0081 | 414668.9933 
98 2132 53.54 906.60 0.059 145.287 1.000 525.272 2096156.0088 | 414810.9941 
151 3397 15.96 423.66 0.038 0.071 11.119 2097854.0075 | 427506.9935 
351 7749 8.97 298.70 0.030 829.408 3.750 329.934 2069017.0080 | 432938.9936 
470 10489 133.32 4 373.31 0.030 177.529 1.400 1 831.073 2090010.0087 | 414803.9927 


y arroyo Grande aguas abajo de San Pedro 
Tlanixco; Granadas, Zacapalco y Tlahuichia 
al sur de Teacalco; y San Jerónimo al sur de 
Zumpahuacán. 


Análisis de costos 
En este estudio se encontraron 534 sitios con 


potencia bruta menor de 250 kW, 41 sitios 
con potencia entre 250 y 1 000 kW, y tres con 


potencia mayor a 1 000 kW. El costo (pesos 
mexicanos) más elevado fue de $43,781,104.00 
para el sitio con una potencia bruta de 1 831.1 
kW, con caída de 133.32 metros y caudal de 1.4 
m?*/s, y el mínimo de $86,357.00, correspon- 
diente a un sitio con caída bruta de 7.36 metros 
y potencia bruta de 1.52 kW. Se determinó 
que el costo total promedio de una PCH es de 
$3,672,636.00 y el costo por kW instalado es 
$52,865.00 en promedio. 
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Figura 10. Distribución de los sitios con potencial micro y mini hidroeléctrico. 


Factibilidad de una PCH al hilo de agua 


Como resultado del análisis de factibilidad se 
obtuvo que el diseño usando el Q,, es el más 
conveniente y se obtendría mayor beneficio por 
la inversión de acuerdo con la TIR, pero el valor 
presente neto indica que el proyecto que aporta 
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mayores ganancias es el obtenido al realizar el 
proyecto con un caudal de diseño igual a Q,, 
como se puede observar en la figura 11. Sin em- 
bargo, dado que se busca que la planta funcione 
la mayor parte del tiempo posible, se considera 
que el Q,, es apropiado para el diseño de plantas 
al hilo del agua. 
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Figura 11. Rentabilidad del proyecto en el sitio 351 con código 7 749, cuenca Chontalcuatlán. 
Conclusiones El método propuesto se puede aplicar 


La planeación global en el sector energético 
facilita el desarrollo de generación de energía, 
considerando fuentes renovables, como las 
pequeñas centrales hidroeléctricas al hilo de 
agua. La metodología presentada en este estu- 
dio pretende dar bases para la futura evaluación 
energética. 

En específico, el modelo hidrológico Soil 
Moisture Accounting es una herramienta re- 
comendable para conocer el caudal medio 
en cuencas no aforadas, siempre y cuando se 
tenga información confiable para realizar la 
calibración y validación de los parámetros que 
se incluyen en el modelo. 

El uso del modelo requiere de cierto cuidado, 
pues pequeñas variaciones en los parámetros 
más sensibles pueden dar resultados conside- 
rablemente diferentes. 

En la cuenca del río Alto Amacuzac se 
encontraron 578 sitios con un potencial bruto 
disponible de 49 227 kW. Esta cantidad repre- 
senta un 33.5% de la capacidad autorizada por 
la Comisión Reguladora de Energía para auto- 
abastecimiento en 2011. Es importante resaltar 
que esta cantidad puede variar dependiendo 
de los criterios considerados para elegir un 
sitio apropiado para la generación de energía 
eléctrica. 


teniendo conocimientos básicos de sistemas 
de información geográfica. No es necesario 
aprender el manejo de un software complejo para 
poder localizar los sitios factibles para construir 
pequeñas centrales hidroeléctricas. 

En la evaluación económica realizada para 
los sitios disponibles, el costo aproximado de 
inversión (en pesos mexicanos) en una PCH 
es de $3,672,636.00 en promedio, y los costos 
estimados variaron de $86,357.00 para el caso 
de la PCH con menor potencial a $43,781,104.00 
para la del potencial mayor encontrado. 

Del estudio de factibilidad que se realizó 
para un caso específico se puede deducir que el 
Q,, es apropiado para el diseño de PCHs al hilo 
de agua de acuerdo con los resultados obtenidos 
de la tasa interna de retorno (TIR). 
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